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ABSTRACT 

We investigate the formation of dust grains in the ejecta of Type la supernovae (SNe la), adopting 
the carbon-deflagration W7 model. In the calculations of dust formation, we apply the nucleation 
and grain growth theory and consider the two cases with and without formation of CO and SiO 
molecules. The results of the calculations show that for the sticking probability of aj — 1, C, silicate, 
Si, and FeS grains can condense at early times of '^100-300 days after the explosion, whereas Fe 
and SiC grains cannot form substantially. Due to the low gas density in SNe la with no H-envelope, 
the average radii of the newly formed grains are generally below 0.01 /Ltm, being much smaller than 
those in Type II-P SNe. This supports our previous conclusion that the radius of dust formed in 
the ejecta is smaller in SNe with less massive envelopes. The total dust mass ranges from 3 x 10~^ 
Mq to 0.2 Mq for aj = 0.1-1, depending on whether or not CO and SiO molecules are formed. We 
also estimate the optical depths and thermal emission by the newly formed dust and compare to the 
relevant observations of SNe la. We find that the formation of C grains in SNe la is suppressed to 
be consistent with the observational constraints. This implies that energetic photons and electrons 
heavily depress the formation efhciency of C grains or that the outermost C-0 layer of SNe la is 
almost fully burned. Finally, we perform the calculations of dust destruction in the SN remnants and 
find that dust grains formed in the ejecta of SNe la are almost completely destroyed in the shocked 
gas before being injected into the interstellar medium. This indicates that SNe la are unlikely to be 
the major sources of interstellar dust. 

Subject headings: dust, extinction — galaxies: abundances — infrared: stars — ISM: supernova 
remnants — supernovae: general — white dwarfs 
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1. INTRODUCTION 

Type la supernovae (SNe la) are considered to be 
thermonuclear explosions of carbon-oxygen white dwarfs 
(WDs). In the single degenerate scenario, the progen- 
itor of an exploding WD grows up to a mass close to 
the Chandrasekhar limit through mass transfer from the 
binary companion and finally undergoes carbon ignition 
near the center (e.g., Nomoto 1982). However, it has 
been a matter of controversy whether the nuclear burn- 
ing front propagates as a subsonic deflagration wave 
(Nomoto et al. 1976, 1984b) or as a supersonic deto- 
nation wave growing from the deflagration (Khokhlov 
1991a, 1991b). Although there are some differences in 
the resulting elemental distributions, both the propaga- 
tion modes can synthesize a significant amount of iron- 
peak elements as well as intermediate-mass elements such 
as Si, S, and Ca (Hoflich et al. 1995; Iwamoto et al. 1999). 
Thus, SNe la are major sources of heavy elements and 
play a critical role in the chemical enrichment history of 
the universe (e.g., Nomoto et al. 1984a; Timmes et al. 
1995; Kobayashi & Nomoto 2009). 

It has been also supposed that SNe la can be possible 
producers of dust grains, especially Fe grains (Tielens 
1998) and would have important consequences for the 
evolution and inventory of interstellar dust (Dwek 1998). 
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This presumption is mainly based on their metal-rich 
compositions of the cooling gas in the ejecta similar to 
those in core-collapse supernovae (CCSNe) for which sev- 
eral pieces of observational evidence for dust formation 
have been reported (see Kozasa et al. 2009 for review). 
In addition, the formation of dust in SNe la has been 
suggested from the analysis of presolar grains extracted 
from meteorites; Clayton et al. (1997) pointed out that 
the isotopic signatures of type X silicon carbide (SiC) 
particles can be explained if SiC grains condense out of 
the explosive helium burning layer in SNe la. However, 
the composition, size, and amount of dust grains formed 
in the ejecta of SNe la and injected into the interstellar 
medium (ISM) have not been fully explored to date. 

It should be mentioned that systematic studies of SN la 
rates have proposed the presence of a prompt component 
of SNe la exploding in a timescale as short as 0.1 Gyr 
after the steUar birth (Mannucci et al. 2005, 2006; DeUa 
Valle et al. 2005; Scannapieco & Bildsten 2005). The in- 
creasing testimonies of such short-lived SNe la have been 
thereafter given by many theoretical and observational 
works (Totani et al. 2008; Aubourg et al. 2008; Hachisu 
et al. 2008a, 2008b; Matteucci et al. 2009; Brandt et al. 
2010, Maoz et al. 2010). Although in the single degener- 
ate scenario SNe la are not expected to generate in metal- 
poor environments as low as [Fe/H] < — 1.1 (Kobayashi 
et al. 1998), prompt SNe la might have occur in the early 
universe at redshifts higher than z = 4, where most of 
observed quasar systems have already reached solar- and 
super-solar metallicities (Juarez et al. 2009). Hence, if 
prompt SNe la could actually produce dust grains, they 
would have made a large contribution to dust budgets 
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not only in our Galaxy but also in high-rcdshift dusty 
galaxies (e.g., Valiante et al. 2009). 

On the other hand, no observation of normal SNe la 
has reported convincing evidence for the ongoing forma- 
tion of dust in the ejecta as indicated by an increase of in- 
frared (IR) continuum, a rapid decline of optical luminos- 
ity, and blueshifts of atomic line emissions; Gerardy et al. 
(2007) observed two SNe la, SN 2003hv and SN 2005df, 
at t ^ 100-400 days since the explosion with the Spitzer 
Space Telescope. However, the observed mid-IR spectral 
energy distributions (SEDs) did not show any evolution, 
allowing the authors to conclude that the mid-IR fluxes 
are dominated by the strong atomic line emission rather 
than thermal emission from newly formed dust. Maeda 
et al. (2009) discovered an accelerated fading of the op- 
tical light curve for the peculiar SN la 2006gz, but it 
remains unclear whether this is responsible for dust for- 
mation. Furthermore, no IR emission from dust formed 
in the ejecta has been seen in the supernova remnant 
(SNR) of Tycho (Douvion et al. 2001) identified as a 
standard SN la (Krause et al. 2008), although recently 
Ishihara et al. (2010) suggested the possible detection of 
thermal emission from the ejecta-origin dust, based on 
the IR observations with the AKARI. 

The fact that the existence of newly formed dust has 
not been distinctly confirmed so far implies that forma- 
tion of dust in SNe la is likely to be inefficient or impos- 
sible because of some physical conditions in the ejecta 
different from CCSNe; SNe la have much higher expan- 
sion velocities of the ejecta and produce much more ra- 
dioactive element ^^Ni than typical CCSNe. A higher 
expansion velocity leads to a lower gas density in the 
ejecta, which prevents dust grains from growing to radii 
larger than 0.01 /im or even makes the condensation of 
dust unsuccessful (Nozawa et al. 2008, 2010). In addi- 
tion, energetic photons and electrons generated from the 
decay of radioactive elements can inhibit the formation 
of dust grains (Nozawa et al. 2008) and also can disturb 
the formation of SiO and CO molecules that affects the 
dust formation processes (Hoflich et al. 1995). Thus, it is 
essential to examine how the formation process of dust in 
the ejecta depends on the type of SNe, by taking account 
of the time evolution of gas temperature and density, the 
formation of CO and SiO molecules, and the impacts of 
high-energy photons and electrons. 

In this paper, we investigate the formation of dust in 
the expanding ejecta of SNe la, with the aim of revealing 
the roles of SNe la in the dust enrichment in the uni- 
verse. In Section 2, we describe the models of SNe la 
and the method of calculation of dust formation. In the 
calculations, wc adopt as a model of SNe la the carbon- 
deflagration W7 model (Nomoto et al. 1984b; Thiele- 
mann et al. 1986), and apply the nucleation and grain 
growth theory to estimate the composition, size, and 
mass of dust that can condense in the ejecta. The re- 
sults of the calculations are presented in Section 3 for 
the two cases with and without the formation of CO and 
SiO molecules. In Section 4, we calculate the optical 
depths and thermal emissions by the newly formed dust 
and apply them to the observational constraints. Then 
we discuss the formation process of dust in SNe la and 
the elemental composition in the outermost layer of SNe 
la. In Section 5, wc calculate the destruction of the newly 
formed grains by the reverse shock and examine whether 
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Fig. 1. — Relative number abundances of elements taking part 
in dust formation in the ejecta of the SN la (W7 model) at 300 
days after the explosion as a function of enclosed mass. The decay 
of radioactive elements is talsen into account, and the abundance of 
each element is summed up over the isotopes. In the calculations 
of dust formation, the original onion-like composition is assumed 
to be preserved with no mixing of elements. 

they can survive to be injected from SNe la into the ISM. 
The summary is devoted in Section 6. 

2. CALCULATION OF DUST FORMATION IN THE EJECTA 

OF SNe la 

2.1. Model of SNe la 

The composition, size, and number density of dust 
formed in SNe depends on the elemental composition in 
the ejecta as well as the time evolution of the density and 
temperature of the gas (Kozasaet al. 1989, 1991; Nozawa 
et al. 2003). In the calculations of dust formation in SNe 
la, we adopt the nucleosynthesis and hydrodynamic re- 
sults of the one-dimension carbon-deflagration W7 model 
(Nomoto et al. 1984b; Thielemann et al. 1986). The W7 
model can account for many of the observed characteris- 
tics of normal SNe la, so it has been extensively used as 
one of the most standard SN la models. In the W7 model, 
the ejecta mass is Mgje = 1.38 M©, and the kinetic en- 
ergy of the explosion is E'kin/lO^^erg = E^^i = 1.3. The 
mass of synthesized ^"^Ni is M{^^m) = 0.6 A/q, which 
is one order of magnitude larger than those (Af(^^Ni) 
~ 0.06 Mq) produced in typical CCSNe. 

Figure 1 shows the number abundances of elements 
relevant to dust formation in the ejecta of the W7 model 
at 300 days after the explosion as a function of enclosed 
mass Mr. In Type II-P SNe (SNe II-P) with massive H- 
envelopes, dust grains can only form in the metal-rich gas 
inside the He core (Kozasa et al. 1989). However, SNe la 
do not have H-rich envelopes, and dust formation would 
be possible in the entire region of the ejecta. We divide 
the ejecta of the W7 model into four layers, according to 
the elemental composition of interest to dust formation: 
the Fe-Ni layer (Mr = 0-0.91 Mq), the Si-S layer {Mr = 
0.91-1.14 Mq), the 0-rich layer {Mr = 1.14-1.28 Mq), 
and the C-0 layer {Mr = 1.28-1.38 Mq). 

In the deflagration W7 model, 0.1 Mq of C and O 
remains unburned in the outermost layer. In W7, the 
initial C+0 white dwarf is assumed to be composed of 
X(i2C) = 0.475, X{^^0) = 0.50, and X{^^Ne) = 0.025 
(where X denotes the mass fraction) with a number ratio 
of C/0 = 1.27. Therefore, we can expect the formation of 
C grains in the outermost C-0 layer, although the C/0 
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ratio in the outermost layer is determined by many cycles 
of recurrent He-shell flashes, thus being highly uncertain. 
On the other hand, in the delayed detonation model, pre- 
existing C atoms can be completely burned, except for 
those in the very outer layers with the expansion veloci- 
ties higher than ~20,000 km s~^ (Iwamoto et al. 1999), 
which prohibits the condensation of large amounts of C 
grains. The formation of C grains and the elemental com- 
position in the outer layer of SNe la will be discussed in 
Section 4.4. 

In the calculations, we consider the onion-like elemen- 
tal composition as shown in Figure 1, assuming that any 
microscopic and macroscopic mixing of elements do not 
occur before the condensation of dust grains. This as- 
sumption is supported by a variety of studies; the spec- 
trum synthesis calculations of a time-series of spectra ob- 
tained during t =5-377 days for SN 2002bo, SN 2004eo, 
and SN 2003du demonstrated that the abundance distri- 
butions in the ejecta of these SNe la are fully stratified 
with some degree of macroscopic mixing (Stchle et al. 
2005; Mazzah et al. 2008; Tanaka et al. 2011). Based on 
the analysis of the emission profiles in the mid-IR spec- 
trum at t ~ 135 days, Gerardy et al. (2007) have also 
shown that the ejecta structure in SN la 2005df is chem- 
ically stratified. Furthermore, X-ray observations of the 
Tycho SNRs indicate that the layered composition has 
been retained in the ejecta even at ^-^450 yr after the ex- 
plosion (Badenes et al. 2006; Hayato et al. 2010 see also 
Koscnko et al. 2010 for Type la SNR 0519-69.0). 

The time evolution of the gas density in the freely ex- 
panding ejecta is calculated as p oc t^^, based on the 
density structure of the W7 model. The time evolution 
of the gas temperature is evaluated by solving the ra- 
diative transfer equation simultaneously with the energy 
equation of the radiation plus gas under the assumption 
of the local thermodynamic equilibrium (LTE, Iwamoto 
et al. 2000). The energy deposition through the decay 
of radioactive isotopes is calculated by giving a constant 
7-ray absorption opacity that can reproduce the behav- 
ior of the bolometric light curves observed for standard 
SNe la. The kinetic energies of positrons are assumed to 
be deposited instantaneously in situ. The energy loss by 
positron escape becomes more important at later times 
when the density in the ejecta is low enough. However, 
at the epochs relevant to this study {t < 300 days), this 
effect can be negligible (Milne et al. 2001). 

Figures 2(a) and 2(b) show the structures of the gas 
temperature and density in the ejecta of the W7 model, 
respectively, at 100 days (thick solid lines) and 300 days 
(thick dashed lines) after the explosion. Figure 2(c) dis- 
plays the velocity distribution of the W7 model. For 
comparison, we also show those in the He core of the SN 
II-P model with an M^je = 17 Mq, an E^i = 1, and 
a massive H-envelope of Mhohv = 13.2 Mq (thin lines; 
Umeda & Nomoto 2002). As can be seen from Figure 
2(c), the expansion velocity of the metal-rich ejecta is 
much higher in the SN la than in the SN II-P, despite 
the almost same explosion energy of both the SNe; in 
the SN la all of the explosion energy is deposited into 
the low-mass ejecta with no H-envelope, whereas in the 
SN II-P most of explosion energy is used to blow off the 
massive H-envelopc. As a result, the gas density in the 
SN la is more than three orders of magnitude lower than 
that in the SN II-P. In addition, the gas temperature in 
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Fig. 2. — (a) Temperature and (b) density distributions of the 
gas in the ejecta of the SN la (W7 model) at day 100 (thick soUd 
lines) and day 300 (thick dashed lines) after the explosion, and (c) 
the velocity structure for the W7 model (thick lines). The enclosed 
mass is normalized to the ejecta mass of the SN la (Mc = Af^je = 
1.38 Mq). For comparison, shown are those for the SN II-P model 
(thin lines) with an Mpr = 20 Mq and an i?5i = 1 (Umeda & 
Nomoto 2002). The mass coordinate for the SN II-P is normalized 
to the mass of the metal-rich ejecta (Mc = Mnocorc — Mmasscut = 
3.34 Mq) within which dust formation is realized (Nozawa et al. 
2003). 

the SN la decreases more quickly than that in the SN 
II-P and drops down to -2000 K at 100 days. This is 
because the lower gas density causes absorption of 7-ray 
in the ejecta to be less efficient and because the thermal 
energy produced at the explosion is quickly lost by the 
adiabatic expansion unlike SNe II-P. 

2.2. Nucleation Rate 

One of the main aims of this paper is to clarify how the 
size and mass of dust formed in the ejecta depend on the 
type of SNe. In order to achieve this goal, we apply a non- 
steady state nucleation and grain growth theory, which 
has been utilized for investigating the dust formation in 
various types of CCSNe in our previous studies (Nozawa 
et al. 2003, 2008, 2010). The theory enables us to es- 
timate the size distribution and mass of newly formed 
dust, given the elemental composition and the time evo- 
lution of density and temperature of the gas. Here we 
present the formula for the nucleation rate, which will be 
useful for discussion in the later section. For comprehen- 
sive description of dust formation calculations, we refer 
the readers to Nozawa et al. (2003). We also note that 
there is another approach to evaluate the amount of dust 
formed in the SN ejecta, which follows the abundances of 
dust precursor molecules by solving a complex chemical 
reaction network (Cherchneff & Dwek 2009, 2010). 

The nucleation rate, taking account of chemical reac- 
tions at condensation, is evaluated by introducing the 
concept of the key species defined as the least abundant 
gaseous species among the reactants. By assuming that 
the key species controls the kinetics of nucleation and 
grain growth (Kozasa & Hasegawa 1987; Hasegawa & 
Kozasa 1988) and taking into account the difference be- 
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TABLE 1 

Chemical reactions for silicate grains considered in this paper 



dust species key species 



chemical reaction 



A/10*K 



B 



(2) 



MgSi03(,) 
Mg2Si04(,) 



Mg(g)/SiO 



model A 

"Sin; 



2U 



MgSi03(,) 
4 Mg2Si04(s) 
-> Mg2Si04(s) 



(3) SiO: 



M. 



SiO 
SiO 



(g) 



:(g) 
(g) 



g(g) -t 

2Mg(g) 
2Mg(g) 



SiO 
SiO 



SiO 



O 



(g) 
(g) 



30 
30 



(g) 



SiO; 



(g) 
(g) 



25.0129 72.0015 

18.6200 52.4336 

37.2400 104.8672 

12.6028 38.1507 



"W" 

(5) 



MgSi03(,) 

Mg2Si04(s) 



MS(g)/«i(g) 



2Mg(g) 
2Mg(g) 



model B 
^i(g) 



- y^(g) - 
+ 40(g) - 
+ 40(g) ■ 

* Sl02(s) 



MgSi03(,) 

■> Mg2Si04(s) 

-> Mg2Si04(s) 



(6) Si02(,) 



Si 
Si, 



(g) 
(g) 
(g) 



Si 



(g) 



20 



(g) 



34.7214 87.7178 

23.4742 60.2918 

46.9484 120.5836 

22.3113 53.8670 



Note. — The key species is defined as the gas species whose coUisional frequency is the least 
among reactants. The subscripts (s) and (g) denote solid condensate and gas species, respectively. 
The Gibbs free energy AG*? for the formation of a condensate from reactants per molecule of the 



key species is approximated by AG^/fcT^^. 



-AjT^ j + _B, where the numerical values A and B 



are derived by least-squares fittings of the thermodynamics data (Chase et al. 1985) in the range of 
temperature of interests. 



tween the temperatures of gas and condensation nuclei 
(Kozasa et al. 1996), the steady-state nucleation rate J| 
for a given grain species j is written as^ 



J- 



T 



cxp 



27(lnS'j)2 



where aj is the probability with which the key species 
stick onto the surface of grains, fij = Attuq ^aj /kTdj is 
the energy barrier for nucleation, k is the Boltzmann con- 
stant, aj is the surface energy, T is the gas temperature, 
Td j is the temperature of the condensation nuclei, and 
and aoj = (3fij/47r)^/^ are, respectively, the volume 
and hypothetical radius of the condensate per molecule 
of the key species. The concentration and mass of the 
key species are represented as ci.j and TOi.j, respectively. 
Note that a factor Ilj , which depends on the partial pres- 
sures of reactant and product gas species except for the 
key species (Chigai et al. 1999), is omitted from the equa- 
tion since we set Ilj = 1 in the calculations (see Nozawa 
et al. 2003). The supersaturation ratio Sj is calculated 

by 



In 5*1 



AG? 



kTr 



E-anP. + iln(^ 



(2) 



where AG° is the Gibbs free energy for condensation of 
the jth dust species from the reactants per molecule of 
the key species, and Pi arc the partial gas pressures of the 
reactants and products. The stoichiometric coefficients 
Vi, which are normalized to the key species, are positive 
for a reactant and negative for a product. 

In the following calculations, we assume Tdj = T, to 
be consistent with a treatment in our earlier studies; dust 
formation in the case where the temperature of the con- 
densation nuclei is not equal to the gas temperature is 
examined in Section 4.2. We take aj = 1 and 0.1 as the 
sticking probability, assuming that it is the same for all 
grain species; for aj < 0.01, dust formation does not oc- 
cur substantially for the SN model considered here. The 

^ There is a typographical error in Eq (3) for nucleation rate in 
Nozawa et al. (2003); cij should be replaced with ^ as is seen 
in Eq (1). 



Gibbs free energy for formation of a grain species is ap- 
proximated as AG^/fcTdj = -A/Td^j + B through the 
least-square fitting of the thermodynamics data (Chase 
et al. 1985). The formation of all possible condensates is 
calculated simultaneously. The grain species considered 
and the data needed for the dust formation calculations 
arc given in Table 2 in Nozawa ct al. (2003). 

2.3. Formation of CO and SiO Molecules 

Since the first detection in SN 1987A (Spyromilio et 
al. 1988; Aitken et al. 1988), formation of CO and SiO 
molecules has been observed in several CCSNe (Spy- 
romilio & Leibundgut 1996; Spyromilio et al. 2001; Ger- 
ardy et al. 2000, 2002; Kotak et al. 2006, 2009; Hunter 
et al. 2009). Formation of these molecules prior to dust 
formation can affect the composition, size, and mass of 
dust grains that condense in the ejecta (e.g., Nozawa et 
al. 2003); carbon (oxygen) atoms bound in CO molecules 
are not available for formation of C-bearing (0-bearing) 
grains, and SiO molecules lock up silicon atoms available 
to form Si and SiC. On the other hand, SiO molecules, 
whose rapid depletion was observed concurrently with 
the onset of dust formation in SN 1987A (Roche et al. 
1991) and SN 2004et (Kotak et al. 2009), are consid- 
ered to be precursors of silicate grains such as MgSiOa, 
Mg2Si04, and Si02 (Kozasa ct al. 1989) as represented 
by the reactions (l)-(3) in Table 1. 

In the astronomical environments hospitable to dust 
formation such as in mass-loss winds of evolved late- 
type stars, it is usually assumed that formation of CO 
and SiO molecules is complete. However, in the ejecta 
of SNe, these molecules are destroyed by impacts with 
energetic electrons and/or charge transfer reactions with 
the ionized inert gas. For SN 1987A, the masses of CO 
and SiO molecules were estimated to be < 5 x 10~^ Mq 
(Liu et al. 1992; Liu & Dalgarno 1995; Gearhart et al. 
1999) and < lO'^ Mq (Liu & Dalgarno 1994, 1996), 
respectively. These masses correspond to formation ef- 
ficiencies less than a few percents, where the formation 
efficiency is defined as the mass fraction of carbon and 
silicon locked in CO and SiO molecules. The disruption 
of CO molecules may enable carbon grains to condense 
out of the oxygen-rich gas in the SN ejecta (Clayton ct 
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al. 1999, 2001; Dcneault et al. 2006; Todini & Ferrara 
2001; Bianchi & Schneider 2007). 

In the ejecta of SNe la that produce much more 
^'^Ni than typical CCSNe, destruction of CO and SiO 
molecules can be more effective clS 0, result of more vig- 
orous radioactivity (Liu 1997). Hoflich et al. (1995) ar- 
gued that CO and SiO molecules can form more or less 
in subluminous SNe la with M(^'^Ni) < 0.3 M© but can- 
not form in normal SNe la with M(^'^Ni) ~ 0.6 Mq. 
Indeed, neither the fundamental nor overtone emissions 
of CO molecule have been detected in ordinary SNe la, 
although Taubenberger et al. (2010) has reported prob- 
able detection of CO emission at ~85 day after the peak 
brightness in superluminous SN la 2009dc. We also note 
that possible formation of SiO molecule is supposed for 
SN la 2005df, but its formation efficiency should be ex- 
tremely low (< 0.01 %, Gerardy et al. 2007). 

In the present calculations, we consider two extreme 
cases for formation of CO and SiO molecules. One is the 
complete formation of CO and SiO molecules, in which 
case the lesser of carbon (silicon) and oxygen is locked 
up in CO (SiO) molecules and SiO molecules act as pre- 
cursors of silicate grains. This is the assumption made 
by a series of our previous works. Although CO and SiO 
molecules can be the dominant coolants in SNe (Liu & 
Dalgarno 1995), we do not take into account the effects 
of cooling by these molecules on the thermal structure of 
the ejecta.^ The other case is no formation of CO and 
SiO molecules, where carbon, silicon, and oxygen take 
part in dust formation as free atoms. It should be men- 
tioned here that even if SiO molecules cannot form in 
the ejecta, silicate grains may condense via the chemical 
reaction paths different from the reactions (l)-(3) given 
in Table 1. Here we calculate the formation of silicates 
by considering the reactions (4)-(6) in Table 1. 

3. RESULTS OF DUST FORMATION CALCULATIONS 

In this section, we present the results obtained from 
the models of calculations outlined in Section 2. In what 
follows, we refer to the cases with and without formation 
of CO and SiO molecules as the model A and model B, 
respectively, with the values of sticking probabilities aj 
attached in the model names. For example, the Al repre- 
sents the model taking complete CO and SiO formation 
with aj = 1. For each of models considered in this pa- 
per, the mass of each dust species formed in the ejecta 
are summarized in Table 2. In the followings, we focus on 
the results in the most optimistic cases of aj = 1 for dust 
formation (i.e., the models Al and Bl), to make a direct 
comparison with our previous studies, where aj = 1 was 
assumed. 

3.1. Cases with Molecular Formation 

Figure 3(a) shows the condensation time tc of dust 
grains formed in the ejecta of SNe la for the model Al, 

^ The cooling by CO and SiO molecules causes rapid decreases 
in the gas temperatures in the C-O layer and the O-rich layer. This 
would lead to the earlier condensation of silicate and C grains and 
their larger average radii than those given in this paper. However, 
the temperature decrease by molecular cooling cannot cause the 
mass of these dust grains to increase appreciably in the case of 
aj =1, where most of the metals available for dust formation are 
incorporated into dust grains (see Section 3.1). We also note that 
if the ejecta were optically thick for CO and SiO emissions, the 
cooling could not operate effectively (Gearhart et al. 1999). 
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Fig. 3. — (a) Condensation times and (b) average radii of dust 
grains formed in the ejecta of the W7 model as a function of en- 
closed mass for the model Al taking complete formation of CO 
and SiO molecules with a, = 1. 



where tc is defined as the time when the nucleation rate 
reaches the maximum (see Nozawa et al. 2003). In the 
outermost C-O layer, C grains are formed at 80-110 days 
from the remaining carbon atoms that are not locked in 
CO molecules. In the O-rich layer, AI2O3 grains first con- 
dense around 100 days, and silicate (Mg2Si04, MgSiOa, 
Si02) grains condense at 120-150 days via the reactions 
involving SiO molecules. Then MgO grains form at 150- 
160 days in the O-rich layer, and FeS and Si grains form 
at 160-250 days in the Si-S layer. Despite the presence 
of the extended Fe-Ni layer, Fe and Ni grains are not 
produced at Mr = 0.15-0.8 Mq; the gas density in this 
region becomes too low for them to nucleate before the 
gas cools down to their condensation temperatures (< 
800 K). Fe and Ni grains can condense at 250-300 days 
only in the innermost region {Mr = 0-0.15 Mq) with the 
highest gas density (sec Figure 2), though their masses 
are very small (see below). 

The above results indicate that a variety of dust grains 
can form in SNe la according to the elemental abun- 
dances in the different layers, given that the stratified 
elemental composition is retained with no mixing. The 
resulting composition of newly formed grains and the or- 
der of their condensation are almost the same as those 
in any types of CCSNe (Nozawa et al. 2003, 2008, 2010). 
This means that the condensation sequence of dust docs 
not depend on the type of SNe, reflecting the fact that 
the elemental composition of the metal-rich gas in the 
ejecta is not greatly different among different types of 
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TABLE 2 

Mass of Each Dust Species Formed in the Ejecta of the SNe Ia 
(W7 model) 



dust species 


Al 




AO.l 




Bl 




BO.l 




C 


5.66 X 10" 


-3 


2.84 X 10 


-4 


3.73 X 10" 


-2 


2.40 X 10' 


-2 


MgO 


3.17 X 10- 


-6 


1.85 X 10 


-9 


9.26 X 10" 


-8 


1.93 X 10" 


-9 


MgSiOs 


7.59 X 10" 


-3 


1.31 X 10 


-6 


1.95 X 10" 


-2 


1.11 X 10' 


-5 


Mg2Si04 


7.01 X 10" 


-3 


1.50 X 10 


-6 


6.08 X 10- 


-3 


6.49 X 10' 


-6 


SiOa 


1.47 X 10" 


-2 


9.94 X 10 


-6 


4.91 X 10" 


-2 


2.21 X 10" 


-3 


AI2O3 


8.18 X 10" 


-7 


7.48 X 10" 


-10 


8.53 X 10" 


-6 


7.71 X 10" 


■10 


FeS 


1.78 X 10" 


-2 


1.53 X 10 


-5 


1.78 X 10" 


-2 


1.53 X 10" 


-5 


Si 


6.30 X 10" 


-2 


3.15 X 10 


-5 


6.40 X 10" 


-2 


3.21 X 10' 


-5 


Fe 


9.52 X 10" 


-5 


1.09 X 10 


-8 


9.52 X 10" 


-5 


1.09 X 10' 


-8 


Ni 


1.48 X 10" 


-6 


2.22 X 10" 


-10 


1.48 X 10" 


-6 


2.22 X 10- 


■10 


Total 


1.16 X IQ- 


-1 


3.44 X 10 


-4 


1.94 X 10" 


-1 


2.63 X 10" 


•2 


Note. — 


The dust mass 


is given in 


units of Mq. The models Al and 



AO.l are the cases where formation of CO and SiO molecules is assumed to 
be complete, with the sticking probabilities of aj = 1 and 0.1, respectively. 



In the model Bl with «, 



1 and the model BO.l with o, = 0.1, it is 



assumed that any molecule never forms in the ejecta. 

3.2. Cases with No Molecular Formation 



SNe. However, the condensation times of dust in SNe la 
are much earher {tc — 100-300 days) than those (<c ^ 
300 days) in SNe II-P (Nozawa et al. 2003). Such early 
condensation of dust is caused by the rapid decrease of 
the gas temperature resulting from the low gas density 
in the ejecta (see Figure 2). 

The low gas density in SNe la also leads to an im- 
portant consequence for the size of dust formed in the 
ejecta. Figure 3(b) exhibits the average radii Oavo of 
newly formed grains as a function of enclosed mass. We 
can see that the average radii of all dust species except 
Fe and Ni grains are under 0.01 fim, being considerably 
smaller than those (cavo ^ 0.01 iim) in SNe II-P (Nozawa 
et al. 2003). As demonstrated by Nozawa et al. (2010), 
the condensation of dust grains in less dense gas demands 
a larger supercooling, and the resulting formation and 
growth of a huge number of condensation nuclei makes 
the average radius of grains small due to the conservation 
of mass of materials available for dust formation. On the 
other hand, in the gas with an even much lower density, 
formation of seed nuclei itself is depressed significantly. 
In this case, the final radius is determined by the compe- 
tition between the collision timescale of the key species 
onto the grain surfaces and the expansion timescale of 
the gas. This is the case for Fe and Ni grains, which can 
acquire the average radii of Oavc 0.01 /.tm 

It should be emphasized that the masses of Fe and Ni 
grains are quite low (< 10~^ -^0) although their radii 
are relatively large. This indicates that SNe la cannot 
be major sources of Fe and Ni grains, contrary to the 
presupposition made in some studies. In the model Al, 
the total mass of newly formed dust is 0.116 Mq and 
the main dust species in mass are Si (0.063 Mq), FeS 
(0.018 Mq), and SiOs grains (0.015 Mq). The mass of 
silicate (Mg2Si04, MgSiOs, Si02) grains is 0.03 Mq in 
total and occupies 25 % of the total dust mass. C grains 
are produced with 5.7 x 10~^ Mq. In the model AO.l 
with aj = 0.1, the total dust mass is 3.4 x 10"'' Mq 
and is dominated by C grains (see Table 2). The average 
grain radii are a factor of 5-10 smaller in the model AO.l 
than in the model Al. However, the condensation times 
of dust grains for the model AO.l are not significantly 
different from those for the model Al. 



Figures 4(a) and 4(b) show the condensation times and 
the average radii of newly formed grains in the model Bl, 
where no molecular formation and aj = 1 are assumed. 
Since dust formation in the inner Fe-Ni layer {Mr < 
0.8 Mq) is not affected by formation of CO and SiO 
molecules, only shown the results in the region of Air = 
0.8-1.38 Mq are in Figure 4. In the model Bl with no 
CO formation, nucleation and growth of C grains in the 
outermost C-0 layer can advance in more carbon-rich 
gas. As a result, the average radii and the mass of C 
grains increase by a factor of 2-3 and 6.6, respectively, 
compared with those in the model Al. Furthermore, be- 
cause of no formation of SiO molecules, more free silicon 
atoms are available for formation of Si grains at Mr = 
1.10-1.15 Mq, leading to a little higher mass of Si grains 
than in the model Al. 

The difference of the results between the models Al 
and Bl also appears in formation of silicate grains. In 
the model Bl, Mg2Si04, MgSiOs, and Si02 grains can 
condense a few tens days earlier than those in the model 
Al. Moreover, the condensation time of Si02 grains is 
earlier than those of Mg2Si04 and MgSiOs, in contrast to 
the case in the model Al. These differences comes from 
the difference in the Gibbs free energies for their conden- 
sation; the reactions (4)-(6) in Table 1 provide the lower 
Gibbs free energies than the reactions (l)-(3). Since the 
lower Gibbs free energy allows the supersaturation ratio 
to rise above unity more rapidly (see Eq (2)) for a given 
temperature evolution of cooling gas, the earlier forma- 
tion of silicate, especially Si02 grains, is realized in the 
model B. The average radii and the masses are enhanced 
for Si02 and MgSiOs, as a result of their formation in the 
denser gas (at earlier times). However, both the average 
radius and the mass are slightly reduced for Mg2Si04, 
compared with those in the model Al; the earlier con- 
densation of Si02 and MgSiOs significantly reduces Si 
atoms available for the later formation of Mg2Si04. The 
total mass of silicate is 0.075 Mq, which is 2.5 times 
higher than those in the model Al. 

Our calculations demonstrate that SiC grains never 
condense in SNe la, regardless of formation of CO and 
SiO molecules. The main reason for this is that in the 
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Fig. 4. — (a) Condensation times and (b) average radii of dust 
grains formed at Mr = 0.8-1.38 Mq in the W7 model for the model 
131 taking no formation of CO and SiO molecules with aj = 1. 



ejecta of W7 model there is no layer where carbon and 
silicon coexist abundantly. However, even if both the el- 
ements are abundant in some layer, as appeared in CC- 
SNe, theoretical calculations (Nozawa et al. 2003; Cher- 
chneff & Dwek 2009, 2010) do not predict the conden- 
sation of SiC grains with radii more than 0.1 /im as 
discovered in presolar grains. Thus, it remains unex- 
plained whether presolar SiC grains originated from SNe 
la (Clayton et al. 1997) or CCSNe (Amari et al. 1992; 
Nittler et al. 1996, see also Deneault et al. 2003; Deneault 
2009). 

In summary, formation of CO and SiO molecules, as 
well as the sticking probability, has substantial impacts 
on the mass of dust grains formed in SNe la; the total 
dust mass for the model Bl (BO.l) without molecular for- 
mation is 0.194 Mq (0.026 Mq) and is much higher than 
0.116 Mq (3 X 10"'' Mq) for the model Al (AO.l) with 
molecular formation. In particular, the difference in dust 
mass is more than two orders of magnitude between the 
models of AO.l and BO.l for which the total dust masses 
are dominated by C grains. In the model BO.l without 
CO formation, the number density of C atoms available 
for dust formation is about five times higher than in the 
model AO.l. Thus, even if aj = 0.1, condensation of 
C grains can still proceed efficiently, and a considerable 
amount of C grains can form in the model BO.l. 

On the other hand, in all of the models considered here, 
the condensation times of dust are significantly early (< 
300 days), and the average grain radii are small (< 0.01 
/im). These conclusions arc not influenced by CO and 



SiO formation and arise from the low density of the gas 
in the ejecta of SNe la with no H-envelope. Our preced- 
ing studies demonstrated that a tiny size of dust grains 
with flavc ^ 0.01 /im condense in envelope-stripped Type 
lb SN 2006jc (Nozawa et al. 2008) and Type lib SNe 
(Nozawa et al. 2010), while relatively large grains with 
Oavo ^ 0.01 /im form in envelope-retaining SNe II-P 
(Nozawa et al. 2003). Thus, we conclude that the ra- 
dius of dust formed in the ejecta depends on the type of 
SNe, and that smaller dust grains are produced as the 
envelope masses of SNe become lower. 

4. DISCUSSION 

As shown in the last section, dust grains of ~ 3 x lO"**- 
0.2 Mq can form in the ejecta of SNe la for the stick- 
ing probabilities of aj = 0.1-1. The formation of such 
an amount of dust grains may show some observational 
signatures such as a fading of optical luminosity and a 
rising of infrared luminosity. However, these events sug- 
gestive of dust formation have not been confirmed by 
any observation of ordinary SNe la. In this section, we 
calculate the opacity and thermal emission by the newly 
formed dust and apply the results to the observational 
constraints. Then we discuss the composition and mass 
of dust formed in the ejecta as well as the elemental com- 
position of the outermost layer in SNe la. 

4.1. Optical Depth by Newly Formed Dust 

We begin with estimating the optical depths by dust 
formed in the ejecta. The optical depth t\ as measured 
from the outermost radius R of the ejecta is calculated 
as follows, 



(r,0=^TA,,(r,i) 



= Y, / dr' / TO,(a)4';*(a)/,(a,r',t)43) 

where fj{a,r,t) is the size distribution function of dust 
species j per volume at a position r, mj(a) — Aira^pj/S 
is the mass of a grain with a radius a and a bulk density 
Pj, and Kf^j{a) is the mass extinction coefficient. The 
maximum and minimum radii of dust (flmaxj" and aminj", 
respectively), as well as fj{a, r, t), are obtained from the 
dust formation calculations in Section 3. Here we assume 
that dust grains formed at a give position distribute ho- 
mogeneously within the corresponding mesh. The bulk 
density pj of each dust species is taken from Table 2 in 
Nozawa et al. (2006), and the optical constants for cal- 
culating K^'^* are taken from the references summarized 
in Section 4.1 in Nozawa et al. (2010).^ Note that the 
optical depths scale as t\ oc t^^ for the free expansion of 
the ejecta. 

The results of calculations show that the total optical 
depths at V band (A = 0.55 /im) are T0.55 = 163 (279) 

^ The bulk density and the source of optical constants for Ni 
grains have not been given in Nozawa et al. (2006, 2010). We 
adopt pNi = 8.87 g cm"'' as a bulk density of Ni grains, which 
is evaluated by using the hypothetical radius agj = 1.377 A of 
the condensate (Nozawa et al. 2003). For the optical constants, 
we assume that the optical properties of Ni grains are the same as 
those of Fe grains. Since the mass of Ni grains formed in the ejecta 
is negligible, this assumption has no impact on the results in the 
rest of this paper. 
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at 300 days for the model Al (Bl); C, Si, FeS, and Fe 
grains are the major donors of the opacity with tq 55 j = 
22 (137), 78 (79), 14 (14), and 49 (49), respectively. At 
300 days, silicate grains only have ro.55,sii = 0.01 (0.02) 
because of their low extinction efhciencies at A = 0.55 
^m, while they provide large optical depths of tiqj ~ 
5.4 (11.7) at A = 10 /im. Note that even for the model 
AO.l (BO.l) with Uj = 0.1, the V-band optical depths 
by C grains reach ro.55 = 2.5 (93) at 200 days, which 
are large enough to affect the optical light curve of the 
SN. The calculation by Sugerman et al. (2006) suggested 
that the clumpy structure of the ejecta can make the 
opacity by dust smaller by about one order of magnitude 
than the homogeneous distribution with the same dust 
mass (see also Ercolano et al. 2007). However, even if we 
consider the effect of dust clump, the V-band opacities 
calculated here are still too high to be consistent with 
the observations. This indicates that the mass of dust 
formed in SNe la must be smaller than our calculations 
predict. 

SNe la produce approximately ten times more ^^Ni 
than normal CCSNe and thus may have stronger radi- 
ation fields in the ejecta. Moreover, dust formation in 
SNe la is realized within the first one year, during which 
energetic photons and electrons are likely to be abun- 
dant inside the ejecta. In such an environment, small 
clusters that are composed of several to a few tens atoms 
are heated up by the impacts with UV-optical photons 
and fast electrons, and their temperature could be higher 
than the temperature of the gas, depending on their op- 
tical properties. The heating of small clusters can lead 
to detachment of atoms from the clusters and thus can 
inhibit the formation of dust. In the next subsection, we 
evaluate the temperature of condensate nuclei, and ex- 
amine how the formation process of dust can be affected 
by the effect of non-local thermal equilibrium between 
the nucleated dust and the gas. 

4.2. Dust Formation in non-Local Thermal Equilibrium 

Condition 

In the calculations of dust formation in Section 3, we 
assumed that the temperature of condensation nuclei Tdj 
is the same as the temperature of the gas T. However, 
Tdj is not generally equal to T, and the resulting temper- 
ature difference can influence the condensation process of 
dust (Kozasa et al. 1996) . In particular, if T < T^j in the 
cooling gas, the supersaturation ratio (In Sj) increases 
more slowly than the case of the local thermal equilib- 
rium with T = Tdj (see Eq (2)), which retards the for- 
mation of dust. Furthermore, the effects of the enhanced 
detachment of monomers from a heated small cluster can 
be seen in Equation (1) via the factors {T/Td^y^^ and 

« Td]- 

The temperature of a condensate is determined by the 
balance between the heating and the cooling. In the 
ejecta of SNe, the possible heating processes of small 
clusters are absorption of UV-optical radiation and col- 
lisions with energetic electrons, while the cooling of dust 
takes place via thermal emission and collisions with 
cooler gaseous atoms. However, the energy transferred 
through collisions with the gas is found to be much lower 
than that through radiative processes. Thus, we neglect 
the heating and cooling by collisions with electrons and 
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Fig. 5. — Time evolutions of the square root of the ratio of the 
gas temperature to the dust temperature (T/Tii j)^^"^ for the main 
dust species: Fe grains at Mr = 0.001 Mq, Si and FeS grains at 
Mr = 0.98 Mq, MgSiOs, Mg2Si04, and Si02 grains at Mr = 1.18 
Mq, and C grains at Mr = 1.33 Mq. 

gaseous atoms. Also, for simplicity, we do not consider 
the effect of stochastic heating in evaluating the temper- 
ature of small clusters; in the SN ejecta with a relatively 
strong radiation field, the temperature of a stochastically 
heated grain distributes symmetrically around the equi- 
librium temperature (Nozawa et al. 2008), which there- 
fore does not affect significantly the IR emission spectra 
calculated in the next subsection as well. 

By balancing the heating by photon absorption and 
the cooling by thermal radiation of dust, the equilibrium 
temperature of dust Tdj (r, t) at a given time t and at a 
given position r is calculated from the equation 

F(r, t)KF,j{a, T,) = 4aBTl^{r, t)Kp^,{a, Td,j), (4) 

where F{r, t) is the flux, ctb is the Stcfan-Boltzmann con- 
stant, and Kpj{a,T) and Kpj{a,T) are, respectively, 
the flux-mean and the Planck-mean of absorption coeffi- 
cient Qx^j{a)\ here we assume that the optical properties 
of small clusters are the same as those of the small dust 
grains. The flux of the radiation field F{r, t) in the ejecta 
has been obtained from the radiative transfer calculation 
in deriving the time evolution of the gas temperature. In 
the SN model considered here, the flux at day 400 ranges 
from F = 1.7 X 10^ erg s~^ cm~^ to F = 8.2 x 10^ erg s~^ 
cm~^ over most parts of the ejecta. Note that the spec- 
tral energy distribution of the radiation field in the ejecta 
at late nebula phases have not been known very well both 
observationally and theoretically. In the calculations, we 
assume that the radiation field has a blackbody spectrum 
with a temperature — Tbb = 5000 K, not depending 
on position and time. 

In Figure 5, we present the time evolutions of the 
square root of the ratio between the gas temperature 
and the dust temperature for the main grain species. 
Note that the calculated dust temperature is not sen- 
sitive to grain radius for the size range (a < 0.01 /im) 
of dust formed in SNe la. The temperatures of MgSiOa, 
Mg2Si04, and Si02 grains have been lower than the gas 
temperature (T/Td,j > 1) since their formation and are 
as low as Tdj = 90-150 K at 400 days. On the other 
hand, Fe, FeS, and Si grains hold - 600 K, 720 
K, and 1000 K at 400 days, respectively, and have kept 
higher temperatures than the gas [T/Td.j < 1). The 
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temperature of C grains becomes higher or lower than 
the gas temperature, depending on position and time. 

Using the calculated time evolution of dust tempera- 
ture and taking account of the effect of non-local thermal 
equilibrium, we performed the dust formation calcula- 
tions for C, Fe, FeS, and Si grains with T/Td^j < 1. As 
the results of calculations, the condensation of Fe, FeS, 
and Si grains is delayed by 50-350 days, and their con- 
densation times are 300-400 day for Fe and FeS grains, 
and 450-500 days for Si grains. In such late times, the 
gas density is too low for dust grains to form efficiently, so 
the resulting grain masses are heavily reduced: ~5 x 10"^ 
Mq for Fe grains, '^4 x 10"'* Mq for FeS grains, and 
~1 X 10~^ Mq for Si grains in both the models Al and 
Bl. Thus, apart from Fe grains formed in the innermost 
ejecta, FeS and Si grains make a negligible contribution 
to the opacity (tq.ssj- < 0.1) after > 300 days. On the 
other hand, the mass of C grains, whose condensation 
times are 5-20 days later than those for the case of the 
local thermal equilibrium, decreases only by less than 
~10 % for all of the models considered in this paper. 
Namely, too high optical depth (tq.ss^c > 10) endowed 
by C grains cannot be mitigated by the effects of the 
non-local thermal equilibrium. 

It should be mentioned here that, using a kinetic ap- 
proach of nucleation, Lazzati (2008) pointed out that 
the nucleation of small clusters with a temperature lower 
than that of the gas can proceed faster in the absence of 
a strong radiation field. However, to examine the sta- 
bility of small clusters in the SN ejecta with a strong 
radiation field, the nanoscale effect such as thermal fluc- 
tuations and the finite energy of small clusters need to 
be treated properly. Such a computation is too time- 
consuming (Keith & Lazzati 2011) and is beyond the 
scope of this paper. 

4.3. Thermal Emission from Newly Formed Dust 

In this subsection, we assess the thermal emission spec- 
tra from dust formed in the ejecta by employing the equi- 
librium temperature of dust given in the last subsection. 
Then, we compare the results with the observed mid-IR 
SEDs and discuss the mass and composition of newly 
formed dust. The derived IR spectra will be also useful 
as a guide to the future IR observations of SNe la. 

The flux densities of thermal emission from newly 
formed dust are calculated with Equations (2) and (3) 
in Nozawa et al. (2008), by taking account of the self- 
absorption of dust emission. In this calculation, we adopt 
the mass of each dust species resulting from the dust for- 
mation calculations under the non-local thermal equilib- 
rium condition for the models Al and Bl (hereafter the 
models Al-non and Bl-non, respectively). The observa- 
tional data are taken from the photometric results of SN 
at ~200 days and ~400 days by Gerardy et al. (2007), 
which may be the only mid-IR measurements of normal 
SNe la at late times. Gerardy et al. (2007) claimed that 
the observed mid-IR spectra are not due to thermal emis- 
sion from dust but are dominated due to line emissions. 
Here, we regard the observed flux densities as the upper 
limits for thermal emission from dust. We take D = 18.7 
Mpc as a distance to SN 2005df (Hamuy 2003). 

Figures 6(a) and 6(b) present the calculated total mid- 
IR SEDs at day 200 and day 400, respectively; the thin 
solid (dashed) lines are the thermal emission spectra in 
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Fig. 6. — Spectral energy distributions by thermal emission 
from newly formed dust for the models Al-non (solid lines) and 
Bl-non (dashed lines) at (a) 200 days and (b) 400 days after the 
explosion. The thin lines are the thermal emission spectra summed 
up over all of the grain species formed in the ejecta before 200 days 
or 400 days. The thick lines arc the mid-IR spectra obtained by 
assuming the absence of C grains. The observational data are the 
photometric results of SN 2005df by the Spitzer and are taken from 
Gerardy ct al. (2007). 

the case where the contribution of C grains is included 
for the model Al-non (Bl-non). In both models and in 
both dates, the thermal emission from C grains formed 
in the inner part of the C-0 layer, as well as from sili- 
cate grains in the inner ejecta, is mostly absorbed, and 
the total IR SEDs are dominated by that from C grains 
in the outermost part of the C-0 layer. (As a result, 
the difference of the IR SEDs between the model Al-non 
and Bl-non is small, even if the mass of C grains for the 
model Bl-non is about one order of magnitude larger 
than for the model Al-non.) It should be noted that the 
calculated flux densities are a few orders of magnitude 
higher than the observed values; if C grains had actually 
condensed in SNe la 2005df, the IR observations would 
have detect their continuum emission. This outcome, as 
well as the very high optical depths, conflicts the obser- 
vations, implying that the formation of C grains in SNe 
la should be extremely inefficient or unsuccessful. 

In order to get further insight into the composition 
and mass of dust formed in SNe la, we also calculated 
the mid-IR spectra in the case where the contribution of 
C grains is removed arbitrarily. The results are drawn 
by the thick lines in Figures 6(a) and 6(b). In the mod- 
els Al-non (solid lines) and Bl-non (dashed lines), the 
IR SEDs at 200 day are produced by silicate grains with 
the mass of 0.03 Mq and 0.075 Mq, respectively. At 400 
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days, in addition to the contribution of silicate grains at 
A > 8 /xm, FeS grains of r^A xlQ-'^ Mq contribute to 
the IR fluxes at A < 8 fim. As can be seen from these 
figures, the calculated SEDs do not contradict with the 
observational data at both day 200 and day 400. There- 
fore, from this comparison, we cannot impose any con- 
straints on the composition and mass of grains species 
other than C grains. In other words, ~0. 03-0. 075 Mq of 
silicate grains can be allowed as a mass of dust formed 
in SNc la within the present observational constraints. 

4.4. C grain formation and the Outermost C~0 Layer 
of SNe la 

The discussion in the last subsection indicates that the 
formation of appreciable amounts of C grains in SNe la 
is incompatible with the existing observational results. 
Thus, the condensation of C grains must be suppressed 
much more than the present calculations predict. In our 
model, one of the factors that produces massive C grains 
is relatively high sticking probabilities of aj > 0.1; if we 
take aj < 0.01, neither C grains nor the other grains can 
condense substantially. Such low sticking probabilities 
are suggested for the formation of C grains in the pe- 
culiar Type lb SN 2006jc; Nozawa et al. (2008) showed 
that Qfj < 3 X 10~^ is required for C grains to obtain 
the dust mass needed to reproduce the IR observations 
of SN 2006jc. They argued that in SN 2006jc C grains 
can form at very early times of ~50 days, when energetic 
photons and electrons are expected to be abundant, so 
the low sticking probability may reflect the effective de- 
struction of small clusters by collisions with photons and 
electrons. Given that the formation of C grains in SNe la 
can also occur at an early epoch of ~100 days, it can be 
possible for precursor clusters of C grains to be destroyed 
signiflcantly. 

Another reason for the formation of a large quantity of 
C grains is the presence of massive (~0.05 Mq) unburned 
carbon in the outer layer. In the W7 model adopted in 
this paper, the deflagration wave halts on the way, and 
the original composition of the WDs with a number ratio 
of C/0 ~ 1.27 is retained in the outermost ejecta. How- 
ever, whether the C/0 ratio exceeds unity is not conclu- 
sive, and the composition of the outer layer of accreting 
WDs is unclear because of many uncertain processes such 
as the ^'^C{a,-f)^^0 reaction rate, convective overshoot- 
ing, dredge-up from the core, accretion rate from the 
companion star, and strength of He-shell flashes (e.g., 
Shen & Bildsten 2009; Y. Kamiya et al. in preparation); 
if the C/0 number ratio is less than unity in the surface 
layer of the WD and if a majority of C atoms should be 
locked in CO molecules, C grains may be unable to form 
efficiently in the cooling ejecta. Thus, observations of CO 
molecules in SNe la can provide meaningful clues to the 
formation process of C grains as well as the pre-explosion 
composition of the WDs, although the existence of CO 
molecules has never been reported for normal SNe la to 
date. 

On the other hand, the mass of carbon deduced from 
the observations is much smaller than in the W7 model. 
From the analyses of the spectra before the maximum of 
the light curves (Marion et al. 2006; Tanaka et al. 2008), 
the mass of carbon in normal SNe la is estimated to be 
at most 0.01 Mq. If, as mentioned in these studies, the 
abundance of carbon is 10-100 times less than oxygen. 



the mass of newly formed C grains may be too low for ob- 
servational signatures of dust formation to be detected. 
Therefore, from a point of view of dust formation, we as- 
sert that the composition of the outermost ejecta is very 
carbon-poor. Such a small carbon abundance might re- 
sult from strong He-shell flashes in accreting high-mass 
WDs (Y. Kamiya et al. in preparation) , or be the results 
of the delayed detonation which burns almost all carbon. 

Apart from the outermost layer, the one-dimensional 
delayed detonation models predict the temperature and 
density structures and the elemental composition well 
similar to those in the one-dimensional deflagration mod- 
els (Iwamoto et al. 1999). This implies that in the de- 
layed detonation models the average radii and masses of 
newly formed grains except for C grains cannot signifl- 
cantly be different from those obtained in this paper. The 
remaining amount of unburned carbon in the delayed 
detonation models is smaller (^^0. 005-0. 03 Mq) than 
0.0475 Mq in the deflagration model W7, but it heav- 
ily depends on the still uncertain transition density at 
which the deflagration turns into the detonation. Thus, 
in order to investigate how the mass of newly formed C 
grains can be changed by the underlying explosion mod- 
els, it will be necessary that comprehensive calculations 
of dust formation are performed for the delayed detona- 
tion models with different transition densities. 

Here it would be worth mentioning dust formation in 
the extremely luminous SNe la. Recently, four over- 
luminous SNe la have been discovered and suggested to 
originate from the progenitors with super-Chandrasekhar 
masses: SN 2003fg (Howell et al. 2006), SN 2006gz 
(Hicken et al. 2007), SN 2007if (Scalzo et al. 2010; Yuan 
et al. 2010), and SN 2009dc (Yamanaka et al. 2009; 
Tanaka et al. 2010; Silverman et al. 2011). Interest- 
ingly, these bright SNe la show much stronger C lines 
in their spectra near the maximum brightness than any 
other normal SNe la. This implies that thick unburned 
C-rich layers are left in their outermost ejecta. Since the 
formation of C grains accompanies the presence of the 
C layer as shown in this paper, C grain may be able to 
condense in these super-Chandra SNe la. Indeed, Maeda 
et al. (2009) reported that SN 2006gz after one year is 
quite fainter than expected from the peak luminosity, 
suggesting that such a fading might be caused by dust 
formation. Taubenberger et al. (2010) also discovered an 
enhanced fading of the light curve around 200 days for 
SN 2009dc. 

The super-Chandra SNe la have high ejecta masses 
and the low expansion velocities. Thus, the density of 
the ejecta is higher than those in normal SNe la, which 
may make dust formation more feasible. On the other 
hand, energy deposition from a large amount of ^^Ni (>1 
Mq) produced in super-Chandra SNe la prevents the gas 
temperature from quickly decreasing down to the dust 
condensation temperature (<2000 K), which may make 
the condensation of dust unfeasible. In order to clarify 
how these competing processes affect the condensation 
time and mass of dust, the formation of dust in super- 
Chandra SNe la should be explored both theoretically 
and observationally. 

5. INJECTION OF NEWLY FORMED DUST INTO THE ISM 

The composition, size, and amount of dust grains in- 
jected from SNe into the ISM are determined by the con- 
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flicting processes between the formation of dust in the 
expanding SN ejecta and the subsequent destruction of 
the dust in the hot gas swept up by the reverse and for- 
ward shocks propagating within the SNRs (Nozawa et ah 
2007, 2010; Bianchi & Schneider 2007; Nath et aL 2008; 
Silvia et aL 2010). In this section, in order to reveal the 
mass of dust ejected from SNe la, we investigate the evo- 
lution of dust in the shocked gas inside SNRs, applying 
the results of dust formation calculations. 

The calculations of dust evolution in SNRs are based 
on the model by Nozawa et al. (2007, 2010). In the mod- 
els, the transport and destruction of dust within SNRs 
are pursued simultaneously with the time evolution of 
the gas temperature and density in the spherical sym- 
metry shocks; by treating dust grains as test particles, 
the deceleration by the gas drag and the erosion by sput- 
tering are self-consistently calculated, according to the 
initial radius and the initial position of each dust species. 
The sputtering yield of each dust species is taken from 
Nozawa et al. (2006). 

Note that the destruction efficiency of dust heavily de- 
pends on the ambient gas density as well as the initial 
grain radius. In the present calculations, we assume that 
the ISM around SNe la is uniform, and consider the hy- 
drogen number density of nH,o = 0.01, 0.1, and 1 cm~'^ as 
the gas density in the ISM, whose gas composition is set 
as the solar abundance. As for the model of dust formed 
in the ejecta of SNe la, we adopt the results for the mod- 
els Al and Bl given in Section 3; the adoption of these 
models is intended to estimate the possible maximum 
mass of each dust species supplied into the ISM. Assum- 
ing that the ejecta interact with the ambient medium at 
30 yr after the explosion, we follow the calculations up 
to a few X 10^ yr. 

Figure 7(a) shows the time evolutions of the total mass 
of dust in SNRs for the models Al (solid lines) and Bl 
(dashed line). We can see that for nn.o > 0.1 cm~'^ 
the dust grains formed in the ejecta of SNe la are com- 
pletely (<10-5 Mq) destroyed before 10^ yr. Such ef- 
fective destruction of dust arises from the following two 
reasons. (1) Since SNe la do not have H-envelope, the 
reverse shock can sweep up the dust formation region 
much earlier (<500 yr) than in envelope-retaining SNe 
II-P (>1000 yr). At such early times, the gas density in 
the shocked ejecta is high enough that dust grains are ef- 
ficiently decelerated and eroded due to frequent collisions 
with the gaseous ions. (2) Since the radii of newly formed 
grains are small (^0.01 /im), they are quickly deceler- 
ated by the gas drag soon after encountering the reverse 
shock. Thus, these small grains are finally trapped in the 
postshock gas and are destroyed by thermal sputtering 
without being injected into the ISM. 

On the other hand, for the ISM density as low as 
nHfi ~ 0.01 cm~^, the deceleration and destruction of 
dust in the hot plasma are quite inefficient, and dust of 
~(l-2)xl0~^ Mq survive in the models of Al and Bl; 
C, FeS, and Mg2Si04 grains with smaller average radii 
(flave ^ 30 A) are trapped in the shocked gas to be com- 
pletely destroyed, whereas a fraction of Si, Si02, MgSiOa, 
and Fe grains with larger average radii (oave ^ 30 A) are 
injected into the ISM. However, the total mass of sur- 
viving dust is dominated by Si grains (see Fig. 7(b) for 
the model Al), whose formation has been expected to be 
greatly suppressed by the effect of the non-local thermal 
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Fig. 7. — (a) Time evolutions of the total mass of the newly 
formed dust within the Type la SNRs expanding into the uniform 
ISM with the hydrogen number density of nn o = 0.01, 0.1, and 1 
cm~^. The solid lines are for the model Al, and the dashed lines 
are for model Bl. (b) Time evolutions of the mass of each dust 
species for the model Al and for n^ Q = 0.01 cm"**. 

equilibrium (see Section 4.2). Thus, neglecting the con- 
tribution of Si grains, we regard ^10"'^ M© (^^5 x 10~^ 
Mq) of silicate grains as a more realistic mass estimate 
of surviving dust for tt-h.o = 0-01 cm~^ in the model Al 
(Bl).^ These results indicate that, even if the ambient 
gas density is very low, the amount of dust ejected from 
SNe la into the ISM is considerably small. We also note 
that typical ISM densities around SNe la are estimated 
to be riH.o = 1-5 cm""^ (e.g., Borkowski et al. 2006) and 
that the ambient density lower than nn.o ~ 0.01 cm""^ 
results in too extended SNRs whose sizes are incompati- 
ble with observed ones of nearby Type la SNRs (Badenes 
et al. 2007). 

It should be mentioned that, as discussed in Nozawa 
et al. (2010), if SNRs expand non-spherically, a part of 
the newly formed grains may be able to evade the de- 
struction and to be injected in the directions that do not 
strongly interact with the ambient medium. Recently, 
a variation of late-time nebular spectra of SNe la has 

^ In the model Bl, in addition to silicate grains of ~5 X 10" 
Mq, C grains of 4 X 10"'^ Mq can survive for ?ih,o = 0.01 cm"'^. 
However, as discussed in Section 4, the formation of C grains is 
problematic, judging from the comparison with the observational 
constraints. Therefore, we do not include here the mass of C grains 
as an estimate of the total mass of surviving dust. 
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been successfully explained by taking into account the 
global asymmetry in the innermost ejecta that is pro- 
duced by deflagration developing from the off-center ig- 
nition (Maeda et al. 2010a, 2010b, see also Motohara et 
al. 2006). Aside from the global asymmetry issue, the de- 
flagration flame is also expected to create a small-scale 
mixing structure (e.g., Ropkc et al. 2007). In addition, 
given that SNe la are the end-products of the binary in- 
teraction, aspherical circumstellar medium might have 
been formed around the exploding WDs, which may re- 
sult in the non-spherical evolution of SNRs. However, 
the analyses of X-ray morphology of young SNRs (Lopez 
et al. 2009, 2010) has showed that the remnants of SNe 
la do not seem to be far from spherical symmetry. Al- 
though the effect of degree of non-sphericity and small- 
scale structure on dust survival is uncertain, we conceive 
here that these effects are not enough to make SNe la 
dominant sources of interstellar dust. 

6. SUMMARY 

We investigate the formation of dust in the ejecta 
of SNe la, adopting the carbon-deflagration W7 model. 
In our calculations, we apply the nucleation and grain 
growth theory, to compare with the results for dust for- 
mation in different types of CCSNe in our earlier studies. 
Wc find that for the sticking probability of aj = 1, vari- 
ous grain species can condense in the stratified ejecta of 
SNe la, although Fe grains cannot form appreciably, con- 
trary to the expectation. The composition of dust grains 
formed in SNe la reflects the elemental composition of 
the gas in the ejecta and is basically the same as those 
in any other CCSNe. 

On the other hand, the condensation times of dust in 
SNe la are much earlier {tc = 100-300 days), and the 
average radius of newly formed dust are much smaller 
(oave ^ 0.01 fim) than those in SNe II-P. This is due 
to the low gas density in the ejecta of SNe la that do 
not have H-rich envelopes. We conclude that the radius 
of newly formed dust depends on the type of SNe and 
that smaller grains condense in SNe with less massive 
envelopes. The total mass of dust that can condense 
in the ejecta of SNe la is significantly affected by the 
sticking probability aj and the formation efficiency of 
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CO and SiO molecules, ranging from 3 X 10"" Mq to 0.2 
Mq for a J = 0.1-1. 

Furthermore, we estimate the temperature of small 
clusters and evaluate the effect of the non-local thermal 
equilibrium on the formation process of dust. We find 
that the non-local thermal equilibrium effect can sup- 
press the condensation of FeS, Si, and Fe grains that are 
otherwise the main donors to the optical depths. We also 
calculate the IR emission spectra from the newly formed 
dust, which would be helpful in discussing the possibility 
of dust formation from the future observations of SNe la. 
On the other hand, from the present observational con- 
straints, we notice that the formation of massive C grains 
is suppressed in SNe la. This implies that the nucleation 
of C grains fails due to the destruction by energetic pho- 
tons and electrons, or that the outermost C-0 layer of 
SNe la are burned by the delayed detonation wave, or 
that the pre-supernova C/0 ratio is exceedingly small as 
a result of strong recurrent He-shell flashes. 

Finally, we examine the survival of dust grains formed 
in the ejecta against their destruction in the SNRs. We 
find that, unless the gas density around SNe la is too low 
{n-a,o < 0.01 cm~'^), the newly formed grains are almost 
completely destroyed in the shocked gas in the course of 
their injection into the ISM. Even if the asymmetric effect 
of explosion is taken into account, SNe la is likely to be 
poor producers of interstellar dust. However, SNe la are 
major production factories of heavy elements which will 
be available for the subsequent growth of the pre-existing 
grains through accretion onto their surface in the dense 
molecular clouds. 
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